A ratchet for heat transport between identical reservoirs 
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A one-dimensional periodic array of elastically colliding hard points, with a non-centrosymmetric 
unit cell, connected at its two ends to identical but non-thermal energy reservoirs, is shown to carry 
a sustained unidirectional energy current. 
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A particle subjected to fluctuating forces far from ther- 
mal equilibrium, in a medium with a vectorial asymme- 
try, must drift in a direction determined by the asymme- 
try, even if the fluctuating forces have zero mean. This 
idea can be seen as a consequence of the "Curie Prin- 
ciple" jll - whatever is not ruled out by symmetries is 
permitted and therefore obligatory. Reviews of the many 
realizations of this concept of a "Brownian ratchet" and 
their relevance to the working of motor proteins in living 
cells as well as to novel particle separation methods can 
be found in j2| ^ . Implementations of the ratchet idea 
generally involve randomly forced particles in periodic, 
non-centrosymmetric potentials. A nonequilibrium sta- 
tionary state is maintained either by having the potential 
alternate between two states or by choosing the random 
forcing not to obey a fluctuation-dissipation relation with 
respect to the damping. 

Can these methods for particle transport without 
chemical-potential gradients be extended to generate en- 
ergy currents without temperature gradients? We know, 
of course, that energy cannot flow spontaneously from 
one thermal bath to another with identical thermody- 
namic parameters, no matter how asymmetrical the con- 
ducting medium linking the two baths. The requirements 
of thermal equilibrium and, hence, of time-reversal in- 
variance, mask the structural asymmetry of the conduc- 
tor. We conjecture, therefore, that an energy current 
should flow through a suitably asymmetrical medium if 
the (identical) baths at the ends were non-thermal. 

Another motivation for the work presented in this pa- 
per is the renewed interest in heat conduction and energy 
propagation in momentum-conserving one-dimensional 
systems [|], |j with a temperature gradient imposed by 
connecting the two ends to heat baths at different tem- 
peratures. For a one-dimensional array of point mass 
particles with elastic collisions, if the masses arc all equal, 
the system is integrable and thermal equilibrium is not 
achieved. Even with unequal masses, momentum conser- 
vation causes the heat conductivity to depend singularly 
lH on the length L of the array, as L^^^ The exten- 
sion of these results to the nonequilibrium regime, where 
a nonintegrable array is connected to non-thermal baths 
at the ends, is clearly of fundamental interest. 

We first present a brief summary of our results and 
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FIG. 1: Inverse current flowing from left to right in an ABC 
array, as a function of array length L. Note that the curve 
is sublinear, i.e. the current decays slightly slower than 1/L. 
Particles recoil from the baths with the exponential speed 
distribution Pijv) of Eq. (|l|), or the heavy tailed distribution 
P2{v) of Eq. (S). The y-axis for the Pi{v) and P2{v) are to 
the left and right respectively. 



then discuss our studies in more detail. We study numer- 
ically the dynamics of one-dimensional arrays of particles 
labelled i = 1, ...L from left to right, with masses {m^} ar- 
ranged in a left-right asymmetric sequence. The particles 
are hard points which interact only upon contact, under- 
going elastic collisions. To the two ends of these arrays we 
attach identical, non-t/iermaZ reservoirs of kinetic energy: 
when the first or last particle collides with the reservoir 
beside it, it recoils with a velocity drawn from a distri- 
bution which is not of the form P{v) oc ?;exp[— /3mw^]. 
Regardless of the type of non-thermal reservoir used and 
the nature of the asymmetry of the array, we always find 
a net macroscopic energy current, confirming our conjee- 
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ture above and demonstrating the robust nature of this 
remarkable "energy ratchet". For one class of arrays, 
the steady-state current J decays slowly with the system 
length L, roughly |^] as the behavior expected of 

a conventional heat current in response to a fixed tem- 
perature difference across a length scale L. For another, 
rather different class of arrays, J is found to be indepen- 
dent of L. In either case, it is clear that the current is not 
merely an edge effect present only near the nonequilib- 
rium reservoirs. Systematics of the dependence on size, 
sequence, and termination can be found in Figs. |^ to ^. 

We now present our study in more detail. We illus- 
trate our general point with two types of arrays: (a) 
the 'ABC array, consisting of three species of particle 
A, B, and C with distinct masses tjia < tub < rric, ar- 
ranged in the periodic but non-centrosymmetric sequence 
...ABCABCABC..., and (b) the 'geometric' array where 
nii+i/mi = r = constant . Since the dynamics of each 
particle is purely ballistic between collisions, we adopt 
an "event-driven" approach: given the present momenta 
of all particles, one knows which pairs will next collide 
and when and, hence, the post-collision momenta and 
energies. If the outermost particle at either end collides 
with the reservoir, it is re-emitted into the system with a 
speed V chosen from a distribution P{v). More precisely, 
the rebound speed is first drawn from P{v) and then 
multiplied by 1/V2m where m is the mass of the par- 
ticle, so that particles recoiling from the reservoirs have 
the same typical kinetic energy regardless of their mass. 
We considered two forms for P{v): 



Pi{v) = exp( 
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Both these distributions are easily generated from the 
standard uniform random variable over [0,1). The for- 
mer is somewhat unusual, since if P{v) is taken to be 
the distribution for particles leaking out of a large ho- 
mogeneous reservoir, the velocity distribution inside the 
reservoir must be p{v) ~ P(v)/v which, for (^, is not 
normalizable. However, such an assumption may be too 
restrictive: the velocity distribution of particles in a gran- 
ular flow bumping against the side walls has under certain 
circumstances been found ||] to be ~ exp(— w), and other 
non-equilibrium methods of energy injection at the walls 
leading to Pi{v) could be conceived of. 

We discuss the ABC array first, for which our studies 
used niA = 1, = 2.5, mc = 6.5, with P\(v) [Eq. (^J, 
and vfiA = 1.0, ttib — 2.5, mc = 9.0. with Pi(y) [Eq. 
(||)]. The mass values were selected in a rough attempt to 
maximize the current flowing, at least for small arrays. In 
both cases, lengths L = 52, 100, 202, 400, 802, 1204 were 
studied, chosen so that the particle at the end of the array 
was always an A. For both choices of bath distribution, 
the steady state has a current J flowing from the left 
to the right of the chain and decreasing slowly with L. 
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FIG. 2: Current across ABC arrays, for different chain 
lengths and terminating particles. The sets of points marked 
A, B, and C in the legend correspond respectively to sequences 
starting with A, B and C particles. The velocity distribution 
at the ends is the exponential P\(v) of Eq.(^. 



Whether we use P\{v) or Pi[v), the decay is no faster 
than the 1 jh of "normal" conduction, and is clearly not 
an exponential decay (see Fig. ||). 

Although the current decays slowly with chain length, 
it is sensitive to the termination of the chains. To illus- 
trate this, we study the current J separately for arrays 
starting with A, B and C, as a function of their length. 
For arrays of short lengths that start with A, the current 
is non-monotone with a cycle-of-three periodicity. A sim- 
ilar periodicity is seen when the chain starts with B or 
with C. There is significant variation between the three 
cases. This is summarized in Fig. ^ where the velocity 
distribution used was Pi(v). 

The sensitivity of the current to the termination of the 
chains might seem surprising, in view of our earlier con- 
clusion (Fig. |l|) that the flow of current is not an end 
effect restricted to short chains. A priori, one could have 
proposed various distinct ways in which the energy cur- 
rent J might have depended on the system length L. 
In the first scenario, repeated collisions in the interior 
would dominate for large sizes. The baths at the ends, 
though nonequilibrium, would merely act as crude ther- 
mostats. It would then be tempting to conclude that 
J ~ exp(— L/^), negligible beyond some thermalization 
length ^. At the other extreme, the non-thermal nature 
of the baths could dominate all other effects and give a 
current independent of L. Our results for for the length 
and termination dependence of the current in ABC arrays 
suggests a scenario between these two extremes: interpar- 
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FIG. 3: Average kinetic energy of particles across chains of 
lengths L = 52, 100, 202. The right hand edges of all the 
chains are made to coincide, to show the approximate L- 
independence of the right boundary region. The length of 
the left boundary region is also approximately independent 
of L. The interior region lengthens as L is increased, and the 
kinetic energy varies slowly here. There is a weak cycle-of- 
three periodicity to the kinetic energy, whose cause is not 
clear. The inset plots the velocity distribution, ln[P(i;)/P(0)] 
vs v'^ , for the particle beside the center. There is significant 
deviation from a gaussian for L — 52, but not for L = 802. 



tide collisions might lead to a rapid thermalization as one 
proceeds from the baths into the interior of the chain, but 
with different temperatures near the two ends because 
of the asymmetric chain and nonequilibrium baths. If 
this scenario is correct, the difference between the ef- 
fective temperatures near the two ends would depend on 
the terminations, but with the terminations specified, the 
energy current in a large chain would be due to the ef- 
fective temperature difference, and would decay with L 
as slowly as for thermal conduction. We have measured 
the average kinetic energy across the chain for different 
values of L. As shown in Fig. |^, we find a rough separa- 
tion into two boundary regions whose size does not vary 
significantly with L, and an interior region (whose length 
increases with L) with a slowly varying average kinetic 
energy. However, surprisingly, the kinetic energy in the 
interior is higher than at the two ends. The velocity dis- 
tribution near the center approaches a gaussian as the 
chain length is increased. 

The existence of a non-zero energy current with ather- 
mal baths can be explained qualitatively. If one consid- 
ers a chain of only two particles, when the heavy particle 
recoils from the bath at its end and collides with the 
light particle, it typically causes the light particle to re- 



FIG. 4: Current in the 'geometric' array, as a function of 
array length L with an end-to-end mass ratio mi/rriL = 5 
(open squares), and as a function of log2(mi/m_t) for L = 100 
(solid circles). There is little variation seen as a function of L 
despite the logarithmic scale, in contrast to the pronounced 
maximum as a function of mi/rriL- The error bars are less 
than the point sizes. The baths emit particles with the long- 
tailed speed distribution P2{v) of Eq. 



bound and itself continues to move forward. The energy 
it was carrying continues to move forward, although dis- 
tributed between both particles. But if the light particle 
carries some energy from the reservoir at its end to the 
heavy particle, most of the energy is reflected back and 
returned to the reservoir it came from. Thus we see that 
energy transport is more efficient in one direction com- 
pared to the other. On the other hand, the light particle 
has much higher typical speeds, and therefore impinges 
against its reservoir much more often. Thus there is a 
competition between the relative inefficiency with which 
the light particle transmits energy to the heavy one, and 
its higher attempt rate. Which one of these wins depends 
on the parameters of the system; we have numerically 
verified (for baths given by Pi{v)) that as the mass ratio 
is changed, the current even changes sign. For (identical) 
thermal baths, the two competing effects cancel exactly. 

We now turn to the second type of array we studied: 
the geometric array, where the ratio of the masses of any 
two subsequent particles is rrii^i/mi — r, where r is inde- 
pendent of i. The geometric array was studied using P2{v) 
only. First, a fixed end-to-end mass ratio rtiL/mi — 0.2 
was chosen, and the chain length varied. (The mass ratio 
r was therefore closer to unity for longer chains.) Chain 
lengths were 52, 100, 202, 400, 802, 1204. A large steady 
state current flows from right to left. As the chain length 
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increases, J asymptotes to a constant value, independent 
of L, in striking contrast to what is seen in the ABC array. 
Since J is (roughly) insensitive to L for fixed mi/mi, 
data was also taken with a fixed chain length of f 00, the 
mass at the right end m^, equal to f , and the mass at the 
left end varied. Not surprisingly, this yields a current 
that is zero when rrii = f , peaks as mi is increased, and 
then decays again as mi is increased further. The peak 
is approximately at mi — 100. 

These results can be justified by recalling that for a 
chain where all the particles have equal mass, energy 
propagation is ballistic. With thermal baths at differ- 
ent temperatures at the two ends, such a chain carries 
an L- independent current. (There would be no current 
with identical baths at both ends for such a chain.) Thus 
for a monotonic chain, where the masses vary only gradu- 
ally, one should measure the length in terms of the "scat- 
tering length" over which transport deviates appreciably 
from the ballistic, i.e. the length over which the mass 
changes by some reasonable factor. If m^/mi is held 
fixed, the "effective length" is then independent of L. A 
chain with a very small effective length would look sym- 
metric, and therefore not carry energy between identical 
reservoirs, while the current in a chain with a very large 
effective length would diminish because of repeated scat- 
tering. Although qualitative, this argument would sug- 
gest that the scattering length should roughly correspond 
to mi/rriL = 100, for which J is maximum. 

The two sets of observations on the geometric chain, 
taken together, lead to a surprising conclusion: if one 
takes a (long) monotonic chain with mi = 10 and 
rriL = 1, and another chain with the same ratio between 
successive masses, with mi = 100 and m^ = 10, the cur- 
rent in the first chain should be sqrt(lO) times the current 
in the second, but when the two chains are connected in 
series (in the correct order), the current will be greater 
than that in either of its two components. 

Several natural extensions of the work presented here 



suggest themselves. For example, consider the case where 
the baths at the ends are in identical steady states which 
are characterized by a parameter A such that A = cor- 
responds to thermal equilibrium. How does the ratchet 
energy current depend on A? Another intriguing issue is 
the effect of inelastic collisions, which are crucial if this 
idea is to be applied to granular matter. An analytical 
understanding of our results would clearly be of great 
interest. Most exciting would be an experimental real- 
ization, perhaps a suitably nano-engineered wire in the 
form of a one-dimensional array with an asymmetric unit 
cell, with the ends driven by cells containing identical, 
exothermic chemical reactions, or perhaps a row of ball- 
bearings with particle masses varying as in our ABC or 
geometric arrays, driven at the ends by gas-fluidization. 
The analysis of the latter would have to include inelastic- 
ity, but it should nonetheless display an energy current. 

To summarize: our work proposes a new mode of en- 
ergy transport, in rough analogy to Brownian ratchets 
but with features quite distinct from those systems. The 
heat transport here requires no manipulation of the bulk 
of the one-dimensional medium, only a structural asym- 
metry. Only the two ends of the "conductor" are held 
in (identical) non-thermal steady states. Neither a time- 
dependent potential nor noise that violates detailed bal- 
ance is required in the bulk. Provided energy is supplied 
at both ends in the same nonequilibrium manner, the 
classical dynamics of the intervening medium extracts en- 
ergy preferentially from one reservoir and sends it to the 
other. No violation of the second law of thermodynam- 
ics occurs: we have simply specified the properties of the 
nonequilibrium reservoirs without saying how they are to 
be maintained that way. 
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We have also measured the kinetic energy and velocity dis- 
tribution across the chains, and find that for fixed m_L/mi 
they depend only on the fractional distance x/L down the 
chain. Unlike for ABC chains, the velocity distribution 
near the center does not approach a gaussian for any L 
and niL /mi that we have studied. 



